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Spherical collapse has turned out to be a successful semi-analytic model to study structure forma- 
tion in different DE models and theories of gravity, but nevertheless, the process of virialization is 
commonly studied on the basis of the virial theorem of classical mechanics. In the present paper, 
a fully generally-relativistic virial theorem based on the Tohnan-Oppenheimer-Volkoff (TOV) so- 
lution for homogeneous, perfect-fluid spheres is constructed for the Einstein-de Sitter and ACDM 
cosmologies. We investigate the accuracy of classical virialization studies on cosmological scales 
and consider virialization from a more fundamental point of view. Throughout, we remain within 
general relativity and the class of FLRW models. The virialization equation is set up and solved 
numerically for the virial radius, j/ v ir, from which the virial overdensity Av is directly obtained. 
Leading order corrections in the post-Newtonian framework are derived and quantified. In addition, 
problems in the application of this formalism to dynamical DE models are pointed out and discussed 
explicitly. We show that, in the weak field limit, the relative contribution of the leading order terms 
of the post-Newtonian expansion are of the order of 10 _3 % and the solution of Wang & Steinhardt 
(1998) (see [1]) is precisely reproduced. Apart from the small corrections, the method could provide 
insight into the process of virialization from a more fundamental point of view. 

PACS numbers: 95.30.Sf, 95.35.+d, 95.36.+X, 04.20.Jb 



I. INTRODUCTION 

The question of how structures form in the Universe is a 
long-standing topic in theoretical cosmology and provides a 
lot of room for discussion. Since the fully non-linear regime 
cannot be accessed analytically, huge N-body simulations 
have been set up to describe structure formation by gravi- 
tationally interacting particles in an expanding background. 
However, these attempts are computationally costly and 
therefore perturbative approaches have been developed in 
order to keep the continuous character of general relativity 
and the FLRW model and make use of methods from fluid 
mechanics. A very simple semi-analytic model of this kind 
is the Spherical Collapse. A spherical, overdense patch 
evolves with the background expanding universe, slows down 
due to its self-gravity, turns-around and collapses. The object 
is stabilized by virialization which prevents it from collapsing 
into a singularity. Despite its simplicity and idealizations, 
this model gives a first insight into the formation of spherical 
halos at all mass scales. The underlying formalism dates back 
to Gunn and Gott in 1972 (see [2]), but has been rediscov- 
ered and continuously extended in recent years (see [1, 3-14]). 

In this work, we are going to use the results of Pace et al. 
(2010) (see [12]) to investigate the process of virialization and 
try to find answers to some remaining questions in this con- 
text. The virial theorem provides a powerful tool to study 
systems in equilibrium, but in order to clarify its role in the 
framework of general relativity, a relativistic version is needed. 
After giving an overview of the classical concepts and the re- 
quirements of relativistic calculations in Sects. II and III, we 
will derive a relativistic version of the virial theorem based on 
the Tolman-Oppenheimer-Volkhoff equation (see [15, 16] for 
the original references) in an Einstein-de Sitter and ACDM 
universe (see Sects. IV, V). In the following, this will be ap- 
plied to the virialization equation in the spherical collapse 
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model and a post-Newtonian expansion will be performed (see 
Sects. VII, VIII). The relativistically corrected results for the 
virial radius and virial overdensity will be discussed and lead- 
ing order corrections are worked out in particular. (Sect. X). 
We will also dedicate a section to the problems occurring when 
this formalism is applied to general DE models and point out 
possible ideas to solve them (see Sect. IX). Throughout the 
paper, we will make use of natural units, i. e. c = 1. 



II. VIRIALIZATION IN THE CLASSICAL 
SPHERICAL COLLAPSE 

In the classical treatment of virialization, there are two major 
ingredients that have to be well-understood. First of all, 
structure formation in the present universe is highly non- 
linear on scales less than 10 Mpc and an evolution equation 
for the spherical patch is needed that takes this non-linearity 
into account. Secondly, the virial theorem has to be combined 
with energy conservation to a virialization condition that 
allows determining the time when collapse stops and the 
system reaches an equilibrium. The key quantities, that are 
assigned to it, are the virial radius normalized to the turn- 
around one, j/vir = Rvir/Rta., and the virial overdensity with 
respect to the background, = p(R v i r ) / pb(«vir) • These are 
general functions of redshift and provide a characterization 
of the equilibrium state of the halo. 

The non-linear evolution equation of a spherical overdensity 
of pressureless dark matter has already been treated in detail 
by many authors (see for example [3, 12, 17]). The resulting 
equation 

describes the non-linear evolution of a spherical, top-hat den- 
sity contrast 8(a) = £ -^ k with respect to the background 
dark matter density pb- E(a) contains all the dynamics of the 
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background cosmological model and is related to the Hubble 
function via the expression H(a) = HoE(a). It will be impor- 
tant to express the virial overdensity A v as a function of the 
turn-around one denoted by (: 



A v = 1 + 8(a vil ) = C 

With X vil = , J/vir 

«t a 



3-vir 
2/vir 



(2.2) 



The virial radius, y v i T , is obtained from the virialization equa- 

QJJ 

tion in which the classical virial theorem T = — -r^- is com- 

2 dR 

bined with the assumption of energy conservation during col- 
lapse. 1 It should be mentioned that energy conservation is a 
very common assumption in the literature and it is not proven 
whether it can actually be applied. Maor & Lahav (2005), as 
well as Wang (2006), pointed out that a homogeneous DE 
component with ti)/-l, —1/3 clearly violates energy con- 
servation between turn-around and collapse (see [8, 9, 13]) 



U + 



dU 
dR 



(2-3) 



In case of an Einstein-de Sitter universe, one simply obtains 
J/vir = | whereas the corresponding virial overdensity (eval- 
uated at collapse scale factor) is given by Av = 18ty 2 ~ 178 
(see [13]). In case of ACDM and dynamical DE models, two 
major classical results have been proposed in the literature: 

• Wang & Steinhardt (1998) (WS afterwards) (see [l]) 2 : 



2/vir 



2 



2 + m - ivv 



Wang (2006) (PW afterwards) (see [13]): 

l-(l + 3w)f l + r, t 



2/vir = 



2- (l + 3w)f («,=-!) 2 + r, t 



(2.4) 



(2.5) 



with the WS parameters r\ v and rjt, and the equation-of-state- 
parameter w given by 



Vt = 2C 1 



fiA^ta 



f2m(<lta) 
/ \ 3 

r)v=K — 

\Ovi 



^A(avir) 
(t^vir ) 



PA 
PA 



The corresponding virial overdensities become functions of the 
collapse (virial) scale factor a c (<z v ir) and reach the EdS value 
asymptotically for small scale factors [a c < 10 _1 ) correspond- 
ing to the matter dominated era. 3 



III. REQUIREMENTS FOR RELATIVISTIC 
CALCULATIONS 

Relativistic treatment of virialization in the same way as done 
in the classical case causes some trouble, because energy con- 
servation is not global in general relativity. A second problem 
has been addressed by Komar (see [19, 20]) stating that iso- 
lated bodies like a spherical halo can only be described exactly 
in asymptotically fiat spacetimes which is generally not given 
in case of FLRW models. A promising way out of these prob- 
lems is assuming that the scale of the halo is much smaller 
than the typical length scale of the background universe given 
by the Hubble radius. If the Killing vector field of the FLRW 
spacetime is considered with respect to this assumption, its 
time-like component greatly exceeds its spatial components, 
allowing to neglect the latter and at good approximation con- 
sider the Killing vector as time-like. 

Fig. (1) illustrates that neighbouring Killing vectors are 
approximately parallel on the scale of the object, which 
means that the radial component of K is extremely small 
on these scales and thus negligible. A detailed quantitative 
analysis of this issue is given in Appendix A. 

From this argument, we can infer three major conclu- 
sions essential for the following considerations: 

• Since r <C Rh, the solutions can be considered as nearly 
asymptotically flat such that isolated objects can be defined 
in GR and the halo mass M is well-defined in the sense of 
a Komar integral. 

• Approximately, one can define the scale of the object as local 
and introduce a local coordinate frame which allows energy- 
momentum conservation, V„T M " = 0, to hold in that frame. 
Since the virialization equation, T v i r + U v i r = I7 ta , represents 
energy conservation, this condition is essential. 

• The approximately time-like Killing vector field K of Eq. 
(Al) allows static solutions of the field equations within an 
accuracy of ~ 0.3 % (see Eq. (A19)). In addition, the typ- 
ical structure of the FLRW metric, especially the fact that 
a global time coordinate can be introduced and a space-like 
foliation with it, infers the orthogonality of K to the under- 
lying hypersurfaces such that the Frobenius condition 4 is 
naturally fulfilled. Thus, we can insert and compare static 
solutions at turn-around and virial redshift in the virializa- 
tion equation T v i r + U v i T = Ut&- Since the virial theorem 
will be applied to the final state which is static by defini- 
tion (within the timescales we consider) , this approximation 
holds in particular for the turn-around being a critical point, 
but not a static one in the exact treatment. 



IV. TOV EQUATION FOR MATTER IN THE 
PRESENCE OF A COSMOLOGICAL CONSTANT 

Relying on the assumptions of the previous section, we can set 
up a spherically symmetric, static spacetime with the metric: 



1 In the following, we will drop the bars over the time averaged 
quantities and implicitly assume time averaging. 

2 Horellou & Berge (2005) ([18]) have proposed a generalization 
due to dynamical DE models, but in the ACDM model both 
results agree. 

3 It has to be mentioned for completeness that this is only true for 
dynamical DE models that have negligible contribution in the 
matter dominated era. Counterexamples are early DE models 
(see results in [12] and references therein). 



j z 

as = — e 



2a(r) 



dr + e 



2b(r) 



dr 2 + r 2 dD. 2 . 



(4.1) 



In addition, we consider a system of two fluid components 
described by 



4 For u> being the corresponding dual vector to K, the Frobenius 
condition states that ui A dui = which turns out to be equiva- 
lent to K being orthogonal to the spacelike sections spanned by 
suitably chosen spatial coordinates, (see [21]) 
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FIG. 1. Killing vector field on the spacelike hypersurface of the universe compared to an object of much smaller radius 



in which 



„ - + = (pt + PtWu v + prg^ (4.2) 



PT — pm + PA, PT = Pm + PA- 

The equation of state of the cosmological constant fluid cor- 
responds to w = = — 1. 

Energy and momentum are locally conserved for the total sys- 
tem as well as for each component separately which means 



matter pressure to vanish at the boundary p(r = R) = 0. 6 
This leads to 

p(/(^) 1/2 -l)+2p A 
P( r ) = — ; —U ( 4 - n ) 



where 



V V T^ = and V ■ v T^ m) = 0, V V T^ W = 0. 

Projecting the conservation equation of the (clustering) mat- 
ter component onto the space perpendicular to the time di- 
rection leads to the relativistic Euler equation for the matter 
fluid 

/wV„T M " (m) = with h a „ = + u a u^. (4.3) 
Working that out, one finds 5 

(p + p) V„m = - V a p + uVuP- (4.4) 

In case of a static configuration (V u p = 0) and with the help 
of Eq. (4.1), Eq. (4.4) reduces to: 



(4.5) 



P + P 



The field equations for the given metric are 



1 - 2b f 1 26' 



-26 / // til | 12 . 

e a — a b + a + 



8ttG(p + P a) (4.6) 



8nG(p - P a) (4.7) 



8ttG( P -pa). (4.8) 



Combining Eqs. (4.6) and (4.7) with (4.5) leads to 



e-» = 



■ - ^(p + PA)r 2 
4tyG (p + p)(p + 3p - 2 P a) 



—P = „ r 

3 1 



^(p + PA)r2 



(4.9) 



(4.10) 



Eq. (4.10) is the TOV equation for the ACDM-model in 
case of homogeneous densities. For solving it, we assume the 



5 In the following sections, we define p m = p and p m = p. 



Inserting Eq. (4.10) and Eq. (4.11) into the hydrostatic equi- 
librium condition (Eq. (4.5)) and integrating with the bound- 
ary a(R) — l/21n(l — AR 2 ) (Schwarzschild-de Sitter solution) 
gives 



e a = (1 - AR 



2\l/2 



3-(f^) V2 / 



(4-13) 



3-/ 

and thus the full metric inside the sphere can be written as 



ds 2 = - (1 - AR?) 



f l-Ar 2 V /2 t 
\\-AB? ) J 



3-/ 



dt 2 



(4.14) 



dr 2 



1 — At'- 



+ r 2 dQ, 2 



which represents the metric of the interior Schwarzschild-de 
Sitter spacetime. 

The well-known exterior Schwarzschild-de Sitter solution 



ds 2 = 



1 _2GM_A^ ]dt 2 + 



dr 2 



■ 1 _2G_M_A£ (415) 



+ r 2 dQ, 2 



matches continuously with Eq. (4.14) at r = R. In this par- 
ticular case, it has to be mentioned that asymptotic flatness 
can only be reached approximately as discussed in Sect. (III). 
Since the scale of the halo is much smaller than the Hubble 
radius (v/Rh ~ 10~ 3 ) we can still assume the object to be 
nearly isolated. We decided to embed the sphere into the 
Schwarzschild-de Sitter spacetime instead of an FLRW space- 
time, because spacetime around the object can be assumed 



In the derivation of Eq. (2.1), the assumption of pressureless 
dark matter is a crucial argument. Nevertheless, for consistency 
with the TOV equation, we have to allow a pressure profile for 
the interior of the sphere. This issue will be discussed below. 
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to be approximately static as well (due to the approximated 
time-like Killing vector field on these scales). In the ordi- 
nary Tolman-Oppenheimer-Volkhoff solution (see [15]), the 
perfect fluid sphere is embedded into vacuum described by the 
Schwarzschild solution. In order to be consistent with this ap- 
proach, the generalization including a cosmological constant is 
embedded into the Schwarzschild-de Sitter spacetime. Never- 
theless, it will turn out that the virial radius and overdensity 
can be predicted consistently with this approach, although a 
dark matter contribution outside the sphere is neglected (see 
Sect. (X) and the weak field limits in Sects. (VII) and (VIII)). 



V. DERIVATION OF THE RELATIVISTIC 
VIRIAL THEOREM 

The pressure profile in Eq. (4.11) contains the radial 
dependence of the pressure in a sphere consisting of a 
cosmological constant fluid and collapsed dark matter em- 
bedded into a background Schwarzschild-de Sitter spacetime. 
When virialization starts, the system can be approximately 
assumed to be in equilibrium which means that it can really 
be described by Eq. (4.10) 7 . In order to derive a virial 
theorem from that, one can take the first spatial moment 
which should usually lead to the virial theorem after time 
averaging. This means that small fluctuations around the 
equilibrium state are averaged out over time such that only 
time averaged quantities (energy expressions) are left in the 
virial theorem. Since the system is already in equilibrium 
and the TOV equation has no time-dependences, the time 
integral drops naturally and all quantities can be interpreted 
as time-averaged. 

Eq. (4.10) is multiplied with r and integrated (averaged) over 
the spacetime volume element (hence taking the spatial mo- 
ment and time-averaging are performed in one step): 



lim — 



— lim — 



p rn 



. (p + p)(p + 3p- 2p A ) 



which becomes 



— 4-7T lim — 



4nG 



I 3 a+b i . i 

pre dtdr 



Ar 2 



(5.1) 



o Jo 

T r R 



= 47r lim — 



tL l 



4ttG 4 (p + p)(p + 3p-2p A ) (5.2) 
3 ' 1 - Ar 2 



dtdr. 



Since all the quantities in the integral do not depend on time, 
the evaluation of the time integral cancels naturally and, while 
interpreting the given quantities as time-averaged, this be- 
comes 



■ 47T 



r R 

I I 3 a + b 

/ pre 
Jo 



= 4tt f 
Jo 



dr 



4ttG 4 (p + p)(p + 3p-2p A ) a+b 
3 r 1-Ar 2 € 



(5.3) 



dr. 



7 The TOV-equation represents the equation of motion of the sys- 
tem in equilibrium. 



Looking at the LHS of this equation in Euclidean space and 
performing a partial integration, we see that: 



- 4tt / pVdr = [-47rpr 3 ] Q + / 12nr 2 dr (5.4) 
Jo Jo 

= J 3pdV (5.5) 
ee 2T. (5.6) 
In consistency with the Euclidian case, we can propose that 



2T = 2T = -4tt 
Consequently, we define: 



rR 

I I 3 a + b i 

/ pre dr. 
Jo 



(5-7) 



2T = 4tt 



/' 

Jo 



4-kG 4 (p + p)(p + 3p- 2p A ) a + b 



dr. (5.i 



3 1-Ar 2 
Inserting Eq. (4.11) and Eq. (4.1) into Eq. (5.8), we obtain 



2T = ^- {1-ArY 2 
3(3 -f) y > 



!' 

Jo 



(2p + 2 PA ) 

. 2 \ 1/2 
1-AR 2 



(5.9) 



(1 - Ar 2 ) 3 / 2 



1-Ar» V /2 f 
1-AR 2 ) J 



dr. 



This is one version of a fully relativistic virial theorem for 
clustering dark matter in a ACDM-background model. Of 
course, other attempts exist in the literature to derive a rel- 
ativistic virial theorem for several purposes. Chandrasekhar 
([22]) derived a post-Newtonian version of the tensor virial 
theorem by investigating the post-Newtonian hydrodynamic 
equations consistently with Einstein's field equations. Bonaz- 
zola (1973) ([23]) has proposed an integral identity consistent 
with general relativity in an asymptotically flat, stationary 
and axisymmetric spacetime. Vilain (1979) [24] considers a 
scalar generalization of the virial theorem to general relativ- 
ity which is valid for spherically symmetric, asymptotically 
flat spacetimes which has been successfully applied to stabil- 
ity studies of perfect fluid spheres. In addition, Vilain's work 
allows to interprete the result of Bonnazzola (1973) geometri- 
cally in the spherical case. Gourgoulhon & Bonazzola (1994) 
([25]) extended the work of 1973 to any stationary, asymp- 
totically flat spacetime in general. Straumann ([21]) proposes 
a virial expression in case of a spherically symmetric, static 
spacetime based on the Komar integral and asymptotic flat- 
ness. Except ([22]), these remarkable results have in common 
that asymptotical flatness is a crucial assumption to the space- 
time which is necessary in order to define isolated objects in 
the sense of a Komar integral (see [19, 20]). We want to em- 
phasize at this point that, strictly speaking, this condition has 
to be valid in our case as well. However, we make use of the 
fact that an isolated object can be approximately defined in 
the FLRW spacetime by assuming the scale of the halo to be 
much smaller than the corresponding Hubble radius. 



VI. RELATIVISTIC GRAVITATIONAL 
POTENTIAL ENERGY 

The modified TOV solution can also be applied to find a rel- 
ativistic expression for the gravitational potential energy of a 
spherical body. The derivation is inspired by the considera- 
tions of N. Straumann (see [21]), but since it is quite technical, 
we refer to Appendix B and quote here only the final result: 
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T 



r 

Jo 



A 2 1 

4-7rr e- 



-dr 



(6.1) 



s/1 - Ar 2 

U = f^ P {l--^ w y (6.2) 

with e denoting the relativistic kinetic energy density which 
is defined in Eq. (BIO) (see Appendix B). 
In case of small gravitational fields given for an object having 
a radius r which is much larger than its Schwarzschild radius 
(this corresponds to Ar 2 <g 1), we can expand Eq. (6.1) and 
Eq. (6.2) to first order: 



T = 



— M — M + 



5 

2neAR 5 



+ 0(A 2 



U = -'- 



3 GM 



-R 5 + 0(A ) 



5 R 



^M PA R 2 + 0(A 2 ). 
o 



The kinetic energy will reduce to the specially-relativistic re- 
sult if gravitational effects are neglected to zeroth order. The 
potential energy contains the Newtonian self-energy of a ho- 
mogeneous sphere as a leading-order term. Thus, classical 
limits can be reproduced showing that Eq. (6.1) and Eq. 
(6.2) are consistently defined. 



i?ta, can be obtained by using the solution of Eq. (2.1) and 
Eq. (2.2): 



c = 



gCgta) 

Pm(ota) 



3M 



l+£(ata), p(-Rta) = 



Pm(ata) = ^Q il m a ta 
=$> Rta, = Ita ' 



(7.5) 



2GM 



H 2 Q^(l + 8(a^)) 



1/3 



Let us consider the classical limit with respect to two assump- 
tions: 

• The sphere's radius is much larger than its Schwarzschild 
radius R s = AR 3 « R, i.e. AR 2 < 1. 

• The cosmological-constant density is much smaller than the 
dark matter density inside the sphere. Since pA is of the 
order of the critical density and p — Avp cr with Av ~ 
95 — 180 8 , this can be assumed safely in our case. 

Expanding Eq. (7.2) to first order in AR 2 and Pa/ p, and 
simplifying it, we end up with 



2/vir = 



Pta + Pa 
|p vir + 2p A 



(7.6) 



Writing this in terms of the Wang-Steinhardt-parameters rjt 
and r/v 9 , this becomes 



VII. VIRIALIZATION EQUATION 

Assuming that energy conservation still holds during collapse, 
the virialization equation states 

[T + U] vil = U tB . (7.1) 

Let us now insert all derived energy expressions and perform 
a change of variable r — >■ y = r/i? ta . After simplifying the 
result, we end up with 



T ■ — 

vir — 



8n 2 GRL 
3(3-/) 



(1 - A vil y 2 iT Rt) 1/2 (2p vir + 2p A ) 2 



/ 

Jo 



( i-A vir;/ ^4 y/ 2 



o (l-A vil y*Rlf 2 3 _/_ 



l--4virJ/ 2 -R t 2 a A 



-dy 



f 

(7.2) 



- 2r) v y vil + (2 + n t ) j/vir -1 = 0. 
Eq. (7.7) can be solved approximately by 



1 



Vvir 



(7.7) 



(7.8) 



2 + vt 

Thus, Eq. (7.2) reduces to the WS limit under the given 
assumptions. 



VIII. POST-NEWTONIAN EXPANSION OF 
THE VIRIALIZATION EQUATION 

The post-Newtonian expansion of the virialization equation 
can be done by simply performing a Taylor expansion of Eq. 
(7.2), however we choose a more elegant way including the 
equation of motion of the collapsing sphere. 10 

We begin with a non-static, spherically symmetric spacetime 
described by 



(7 v ir = 4irp viT R t . 



a / 

JO 



v 1 1 



E/ta = 47rp ta i? t 3 a / y 2 I 1 



Jo 



with the definitions 



A v ir — 



8tvG 



v/1 - A vil y*R\ 



8nG 



dy (7.3) 



dy (7.4) 



(p V ir + pA) , A ta = -^T— (Pta + Pa) , 



3M 



f - 1 _ r, - 

' Pvi,' ^-4^^ 



Pta = 



3M 



Eq. (7.1) has to be solved numerically for y vlI at different 
redshifts (see Fig. 2 for the results). The turn-around radius, 



.1) 



and use again the energy-momentum tensor of an ideal fluid 
with two components 



rp _ rp(m) , rp(A) 
± li-V 1 fj,iy i 1 fj,iy 

= (Pm + PA + Pm + Pa)u^U v + (p m + Pa)9^- 

The A-component satisfies an equation of state given by 



(8.2) 



8 See Pace et al. (2010) ([12]) for their results in the ACDM case. 

9 See Wang & Steinhardt (1998) ([1]) and Sect. II 

10 This equation was first derived by Misner and Sharp (1964) (see 
[26]). 
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Pa = -Pa- (8.3) 

Consider a comoving reference frame in which the four velocity 
u has the components 



expanded as well. In our case, spacetime is described by the 
TOV metric given by 



ds 2 = - (1- AR 2 ) 



'3-(^) 1/2 / 



2 . dr 2 



3-/ 



dr + 



1 - Ar 2 



The relativistic Euler equation for the (clustering) matter 
component /i QM V I ,T Ml/ ' (m) = states (in that frame) 11 : 



P 



e a = ^— = -. (8.4) 

p+p p+p y 

If we combine this relation with the field equations for the met- 
ric, we obtain the relativistic equation of motion for a spher- 
ically symmetric object (first derived by Misner and Sharp 
(1964) in [26] and applied by Collins (1978) in [27]) 



d ( dr 



V dt Vdt 



1 dp ( , 2 .2 

— — -r 1 + F 
p + pdr \ 



2GM(r) 
r 



(8.5) 



In the presence of a cosmological constant A with an equation 
of state like Eq. (8.3), this is slightly modified 



d ( dr 



y dt{ y dt 



1 dp I , ■ 

-3- 1 + y r 

y dr 



(p + PA)r 



4ttG 



(8.6) 



(p — 2/9a) r — 4-irGpr. 



In case of equilibrium, Eq. (8.5) and Eq. (8.6) reduce to the 
TOV equations with or without a cosmological constant: 



/ _ 4ttG (p + p) (p + 3g) 
P ~ 3 T 1 - ^fpr 2 

_ ' = 47rG r (P + P) (P + 3 P - 2 P A ) 



1 - 



8ttG 



(p + Pa) r 2 



(8.7) 
(8.8) 



If only small oscillations of the system around its equilibrium 
are considered, we can assume that r 2 /c 2 <C 1. Terms of this 
kind will be neglected in the following. After performing a 
Taylor-expansion up to the first post-Newtonian order O(^) 
and inserting the zeroth-order expansion of the TOV-equation 
(Eq. (8.8)), 



, 47rGr , „ . ^ ( 1 
-P = -^-p(p-2 P a) + 



(8.9) 



we arrive at 



2 .. 

py r = 



4irGr 



P(P~ 2pa) 



-( 



4nGr 



V 3 



(p - 2p A ) + 4irGpr p 



(8.10) 



32tt 2 GV , 2 2 , 
g P (P - PPA - 2p A/ 



Since the derivation of the virial theorem requires an integra- 
tion over the spacetime volume element, the metric has to be 



+ r 2 dQ, 2 



(8.11) 



where 



. 87rG , si.-, 

A =— (P + PA), / = 1 
An expansion up to C(l/c 2 ) leads to 



2p A 
P 



ds 



1 - 2AR 2 



f 



(3-/) 



A(R 2 



r 2 ))dt 2 



(8.12) 



+ (1 + AR 2 )dr 2 + r 2 dQ, 2 



The canonical volume form becomes 



-g AVt 



1 - AR 2 \ 
1-Ar 2 ) 



l^ArtV' 2 f 
1-AR* ) J 



3-/ 



dV T 



(8.13) 



with dVr = dt A dV being the total volume element for a flat 
spacetime in spherical polar coordinates 

dV T =r 2 sin(9-di A dr A d<9 A #. (8.14) 

In the following, we will also apply the definition of the canon- 
ical volume form of the spacelike 3-hypersurface E described 
by the spatial coordinates 

dV s = yj-g\v r 2 sin(9-dr A d0 A #. (8.15) 

Taking the first spatial moment (multiplying with r and inte- 
grating over the spatial volume) leads to the post-Newtonian 
version of Lagrange's identity (see [27]): 

J py 2 rrdVs = - J rp'dVs - J 47r ^* r (p - 2p A ) pdVs 



I 
I 



327r 2 G 2 r 4 , 2 



(p 2 - PPA - 2pi) pdVs- 



(8.16) 



In analogy to the classical case, we interpret 



ld 2 I r _ f 2.. , v 

2T=- J rp'dVs. 



12 I r is defined to be the relativistic generalization of the classical 
moment of inertia. For a homogeneous sphere it is classically 
11 In the following, we define again p m = p, p m = p. defined by I = | J v pr 2 dV. (see [9, 27] 
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Using these definitions, Lagrange's identity becomes the fa- 
miliar expression 



1 d 2 I r 

2~W 



2T- 



4nGr' 



(p-2p A )pdV s 



47i ^ r (p - Pa) + ^TiGpr 2 ) pd\ : 

32iv 2 G 2 r 4 , 2 „ 2\ j„ 
a (P - PPA - 2p A ) pdV s . 



.17) 



Dropping the corrections in 1/c 2 , the classical version of La- 
grange's identity is 



1 d 2 I r 



2T + U m -2U A . 



(8.18) 



2 dt 2 

Performing the time average will lead to the post-Newtonian 
virial theorem, because motions like oscillations around the 
equilibrium configuration are averaged out. Since we have to 
apply 



lim 

T-)-tx> 



A 

Jo 



) e a dt, 



.19) 



the volume element of the averaged form changes dVs = 
e b dV — > e a+b dV = y^gdV while the time integration is per- 
formed. 

Dropping all terms of C(l/c 4 ), the post-Newtonian virial the- 
orem is: 



2T 



AtyGt 2 



(p-2p A )pdV 



- I ( '"J,"' '■ P 2p A ) ' ^Gpr 2 ) (1) 



/ 



9 

2nGr 



+ 1 o (p 2 -PPA-2pi)pdi/ (ii) 

(P - 2p A ) (l + A (r 2 - i? 2 )pdV (III). 

(8.20) 

It can be seen that the correction terms contain 

• pressure contributions (I), since pressure acts as a source 
of gravity 

• backreaction terms (II) between the fluid components 
and the geometry of space (due to the non-linearity of GR) 

• metric expansion terms (III), since a non- vanishing 
energy-momentum tensor changes the metric (due to the 
field equations) 

The potential energy given by Eq. (6.2) can be expanded in 
the same way: 



U = -4np 



^ 4 , 3 .2 6 i i 

— r + -A r ] dr. 



.21) 



Performing the angular integration for the kinetic energy ex- 
pression leads to 



8n 2 Gr 4 



L ~ 3" ' ~ 2/ ' v ' '" ir 



+ 



' ' 8n 2 Gr 4 
3 

„3,-<2„6 



(p — 2pa) + 8ir 2 Gpr 4 ^j pdr 



+ 



f R 647r a G 2 r 6 , 2 2 x , 

J g (P - PPA - 2p A ) pdr 



/ 



, 32tt 3 GV . „ , . , / / 

+" / » (P - 2pa) (p + Pa) p 1 + 



3-/ 



Now we rewrite some variables 13 

r = y ■ i?ta, -Rvir = «/vir ■ -Rta 

and the virialization equation becomes 



[T + U] vil = Uu 



(8.23) 



with the terms 



-47rp v 



Jo 



+ 3 -A 2 viiy e Rt ) Rldy (8.24) 



V + l^tV^a ) ^a*/ (8.25) 



and 



/•Vvir 

Jo 



Gy 4 Rl 
3 

.2/-X..4 r>5 



(pvir — 2pA) pvirdy 



^ (Pvir - 2p A ) + 87T 2 GpvirJ/ 4 ^L ) Pvirrfj/ 



Jo V 3 

^ 647v 3 G 2 y 6 RL , 2 



/ 

Jo 

rV 

Jo 



9 

t 3/-y2„,4 D 7 



Pvir - PvirpA - 2p A ) p vir dy 



32n*G 2 y*RL , 3 2 



(Pvir - Pvir PA - 2p vir p A ) 



1 + 



^3/) (y 2 - 2/vir) ■ 



.26) 



In addition, we have applied the following definitions: 



3M 3M _ 8ttG 

Pvir_ ^y'tirRV Pta ~47rfl t 3 a ' Avir " 3 tPvir + p A ), 
8ttG , 

^4ta = — ^— (Pta + PA) • 

As for the fully relativistic version, Eq. (8.23) has to be solved 
for j/vir numerically. 



IX. THE RELATIVISTIC FORMALISM AND 
DYNAMICAL DE MODELS 

Even though we have spent some effort to generalize our 
method to dynamical DE models, certain problems occur 
which will be described in the following: 

Consider a two component fluid described by 



= (P + p) u^u v + pg^ 
TjS ] = (pq +Pq)u i _ 1 u 1 , +pqP^ 



T _ rp(m) , rp(Q) 
± ILV — ± nu l~ 1 ill/ 



(9.1) 
(9.2) 

(9.3) 



where the densities p and pQ are assumed to be constant 
and the quintessence component has an equation of state 



(r 2 - R 2 ) dr. 



(8.22) 



13 i?ta is again calculated using Eq. (7.5). 



S. Meyer et al. 



Relativistic virialization in the Spherical Collapse model 8 



Pq — wpQ with constant w. Energy-momentum conservation 
is separately fulfilled for each fluid component 

V M T Ml/ ' m = (9.4) 

The static, spherically symmetric field equations for this set- 
up are 

G M „ = 87rGT M „ 

^-e-*(i-^)=ftrG(p + p Q ) (9.6) 

-^+ e "(^ + v) =8 " G(p + Ww) (9 ' 7) 
e~ 2b (a" - a'b' + a' 2 + = 8ttG (p + wp Q ). (9.8) 

It turns out that, in case of w 7^ — 1, Eqs. (9.6) - (9.8) are 
no longer consistently solvable and lead to contradictory solu- 
tions (see Appendix C for a detailed proof). This means that 
a model based on 

• staticity 

• spherical symmetry 

• matter model (two component fluid consisting of cluster- 
ing dark matter and homogeneous DE with equation of state 
PQ = wpq) 

does not lead to a consistent description. Static, spherically 
symmetric solutions with w 7^ — 1 can therefore only be 
given approximately. In order to achieve exact solutions in 
the general case, we need to drop at least one of the model 
assumptions. Since spherical symmetry is dictated by the 
spherical collapse model and the matter model and we want 
to fix the matter model, we can only stick to time-dependent 
problems and drop staticity. 14 

A physical explanation for this constraint on w is local 
energy-momentum conservation. Since we assume that it has 
to be valid for each component separately and the field equa- 
tions are constructed in a way that energy and momentum 
are locally conserved, a static quintessence component with 
w 7^ — 1 violates this requirement. 

Consider a perfect fluid with T M „ = (pq +Pq)m m u„ + PQg M „ 
which has to obey local energy-momentum conservation 

V M T M "' g = 0. (9.9) 

Projecting this onto the 3-space perpendicular to the mean- 
fluid flow, we obtain the relativistic Euler equation 

(<?«„ + u a u v ) V M T'" / ' Q = (9.10) 

which leads to 

(PQ +Pq) V u u = -grad pq - uV u pQ. (9-11) 



In case of a static configuration (u = (1,0,0,0), pci(r,t) = 
PG>{r)) and the metric ansatz from Eq. 4.1, we are left with 

-Pq = a ' (PQ +PQ) (9-12) 
which becomes for the given equation of state pq — wpQ : 

= o'p Q (l + w). (9.13) 

Since a' is non-zero in general, we have to require w = — 1 
in order to satisfy local energy-momentum conservation. 
Thus, a quintessence fluid with constant density and 
violates local energy-momentum conservation in case of a 
static configuration. 

Energy and momentum are certainly conserved for a time- 
dependent DE density which scales like p Q = p^ a~ :i(1+w)15 . 
In this case, the only time-independent DE density is the 
cosmological constant representing the only possible static, 
quintessence fluid configuration with constant equation of 
state that satisfies local energy-momentum conservation. In 
our approach, we restricted ourselves to homogeneous DE that 
evolves independently from the matter component. In case of 
clustering DE (see [5, 8, 9, 11]) or even models considering 
interactions between the two fluids (see [7, 14]), our formal- 
ism will significantly change and might allow an application 
to these models. 



X. RESULTS 

The relativistic virialization equation and its post Newtonian 
expansion are solved for the virial radius. Fig. (2) shows the 
virial radius and overdensity as a function of the collapse red- 
shift z c for three different halo masses in an EdS and ACDM 
cosmology. All quantities at virialization are evaluated at 
the exact virial redshift z vir . It is a common approximation 
in the spherical collapse model to insert the potential energy 
evaluated at the collapse redshift. As already investigated 
analytically by Lee & Ng (2010) (see [6]) the result of the 
virial overdensity changes significantly 16 by inserting the 
exact virial redshift instead. We have developed and applied 
an iterative method to obtain z v i r numerically and the 
results of Lee & Ng are nicely reproduced. The derivation is 
postponed to Appendix D. 

It can be seen clearly in both figures that the Wang- 
Steinhardt-limit is precisely recovered. This is expected, be- 
cause the expressions for the potential energy and the ki- 
netic energy derived in Sect. V contain the WS solution as 
limit to zeroth order. Since the halo mass does no longer 
cancel out naturally on both sides of the virialization equa- 
tion, the spherical collapse becomes mass-dependent. There- 
fore each result is plotted for three different masses, namely 
1O 13 M , 1O 14 M and 10 15 M Q . Nevertheless, it can be seen 
that the results for M = 1O 13 M and M = 1O 15 M differ 
from each other by only 10~ 3 %. 

Since, averaged over sufficiently large timescales, the final 
state of a virialized cluster is static 17 , it can be considered 
as a homogeneous, static perfect-fluid sphere such that the 
Tolman-Oppenheimer-Volkoff solution can be applied. We 



It has to be mentioned for completeness that there exist static ex- 
terior solutions describing a Schwarzschild black hole surrounded 
by a Quintessence fluid (see [28, 29].). However, in these cases 
Pq is constrained to be radial-dependent: Pq(t) ~ r -3(l+iu) 



Of course, this is only true for constant equation-of-state param- 
eter w. In the general case, pq evolves like pq = p>Q g(a) with 
g(a) = cxp (—3 f" (1 + w(a'))dlna') (see for example [4]) 
In the EdS case Ay(z c ) ~ 178 changes into Ay(2 v j r ) ~ 146. 
Oscillations around the virial radius can be expected to average 
out. 
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have constructed the virialization equation based on the TOV 
solution instead of using the classical approach from Fried- 
mann's equations. A few points have to be discussed concern- 
ing this method: 

• As an important approximation, we have assumed that the 
Killing vector field K of the FLRW universe is time-like 
on halo scales. Since the final state is static anyway, this 
is most relevant for the turn-around which is described by 
a static solution as well. As shown in Appendix A, the 
space-like component of K is by at least two orders of mag- 
nitude smaller than the time-like one. Looking at the small 
corrections in first post-Newtonian order in Tab. I being 
five orders of magnitude smaller than the classical term, it 
remains an open question, whether a time-dependent ap- 
proach that does not obey that approximation, would have 
a non-negligible effect on the results. 

• As can be seen in Tab. I, the normalized contributions to 
the first post-Newtonian order are of about 10" 3 % being 
almost independent of the type of contribution. This can 
be expected due to a simple estimate. Let us assume a 
typical massive galaxy cluster with a mass of 1O 15 M and 
a virial radius of f Mpc. The gravitational potential 4>/c 2 
being the ratio of its Schwarzschild radius and virial radius 
has the value 



Einstein de-Sitter: 



GM 

L. l tvir 



5- 10" 



(10.1) 



Thus, the post Newtonian terms which are of the order 
(4>/c 2 ) 2 must have absolute values of about 10^ 10 which 
corresponds to relative contribution of 10~ 5 (10~ 3 %) with 
respect to the classical term. 

The key question remains why the relativistic calculation re- 
duces to the WS limit instead of the result of PW. The ansatz 
of a static, spherically symmetric metric reduces Einstein's 
field equations to a coupled system of three ordinary differ- 
ential equations with respect to the radius r (see Eqs. (4.6) 
- (4.8)). Eq. (4.6) is already decoupled from Eqs. (4.7) and 
(4.8) such that the rr-component of the metric is constrained 
to be 



_ 26 _ 

Qrr — 6 — 



1 



. 8nG . 



!_ Ara . — 3 v.-PAj- (10-2) 

It has to be mentioned that Eq. (10.2) does not contain any 
pressure term. Straumann's self energy expression derived in 
Appendix B is based on the number density n(r) of dark mat- 
ter particles that is integrated over the covariant volume ele- 
ment restricted to the hypersurface £ spanned by the spatial 
coordinates: 



Finally, the result becomes 



VI - Ar 2 



U : 



f R 2 f 
I 4nr p I 1 



Vl - Ar 2 



drd6dcj>. 



dr 



(10.3) 



3 GM _ 4tvG 
~5~ R 5~~ 



(10.4) 



Mp A R 2 + G(A 2 ) 



which perfectly reproduces the potential energy used by 
Wang & Steinhardt in their model. 

The TOV equation for the pressure profile p(r) is a gener- 
alization of the Newtonian hydrostatic equation. Pressure 
is contained as an additional source of gravity, but the 
post-Newtonian expansion of the resulting virial theorem in 
Eq. (8.20) shows clearly that the limit of WS is reproduced to 
zeroth order. Horellou & Berge have shown in 2005 (see [18]) 



pressure term 1.274132394 • 10~ 5 
backreaction term 3.18533352 • 10~ 5 
metric expansion 9.5560225 • 10~ 6 



ACDM: 



pressure term 1.2893437 • 10~ J 
backreaction term 2.7275899168 • 10~ 5 
metric expansion 7.073241606 • 10~ 6 



TABLE I. The three relative post-Newtonian contribution 
terms with respect to the classical Newtonian term are con- 
sidered in EdS and ACDM cosmology for z c — and M = 
1O 15 M . 



that, from a classical point of view, the Wang-Steinhardt- 
solution is exactly valid in the ACDM case. This means that 
for the classical self-energy expression of a matter sphere in 
a ACDM background inserted into the classical virialization 
equation (see Sect. II) the solution of Wang & Steinhardt is 
recovered. 

To conclude this section we can state the following: Based 
on the assumptions of a static, spherically-symmetric space- 
time and separately fulfilled energy-momentum conservation 
equations for matter and cosmological constant 18 , the field 
equations and the resulting TOV equation provide a relativis- 
tic virial theorem that contains the Wang-Steinhardt-result as 
limit to zeroth order. It remains an open question whether a 
non-static approach to relativistic gravitational collapse that 
might also be adapted to dynamical DE models would still 
recover this result. 



XI. REMARKS ON THE REQUIRED 
PRESSURE PROFILE 

As already mentioned by Oppenheimer and Snyder in 1939 
(see [30]), there does not exist any static, spherically- 
symmetric solution of the field equations with vanishing pres- 
sure. If no positive pressure profile is present, nothing will 
prevent the spherical object from collapsing into a singularity. 
Therefore, in the classical spherical collapse, the virialization 
condition is introduced to define an equilibrium. In our case, 
gravity forces the system to build up a pressure profile to pre- 
vent the sphere from collapsing into a singularity. In fact, 
the process of virialization must convert ordered motion from 
the collapse into unordered motion and the kinetic energy as- 
sociated with the unordered motion corresponds to pressure. 
Thus, an equilibrium is reached by a positive pressure profile 
that can be related directly to the mean kinetic energy of the 
system (see Eq. (5.7)). This might be a first step towards a 
more fundamental theory of virialization avoiding of the in- 
terpretation of an "enforced" equilibrium. 
The non-linear density evolution equation (see Eq. (2.1)) 
is still based on extended Newtonian theory. 19 In order to 
achieve full consistency, a relativistic evolution equation for 
either the density or the radius is recommended. Nevertheless, 
these equations will retain their full validity if dark matter is 



18 The cosmological constant fluid trivially fulfills energy momen- 
tum conservation, since the continuity equation reduces to pa = 
0. 

19 Velocities and gravitational potentials are assumed to be small 
compared to c 2 and pressure is included as a source of gravity. 
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Virial radius for an EdS model evaluated at z=z„: 



Virial radius for a LCDM model evaluated at z=z. (i 
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FIG. 2. Virial radius and virial overdensity obtained by the relativistic virialization equation and its post-Newtonian expansion 
as function of collapse redshift for EdS and ACDM cosmologies. The upper panels show the virial radius for three different halo 
masses obtained by solving the full relativistic virialization equation (Eqs. (7.1)-(7.4)) and its post-Newtonian expansion (Eqs. 
(8.23) - (8.26)) using the two cosmological models. The lower panels show the corresponding virial overdensity obtained by Eq. 
(2.2). The classical solutions of Wang & Steinhardt (1998) (cyan line) and Wang (2006) (purple line) are plotted for reference. 
The region close to z c = has been scaled up to sufficiently small redshifts in order to illustrate the dependence on the halo 
mass. 



assumed to be pressureless during the evolution and that the 
pressure profile is built up instantaneously at virialization 20 . 
The appearance of a pressure profile in this context is a di- 
rect consequence of the spherical collapse model itself. Strictly 
speaking, virialization is a tool in the spherical collapse model 
in order to achieve an equilibrium state. The pressure is di- 
rectly related to the mean kinetic energy of the system which 
is again related to the potential energy by the virial theorem. 
Since the last two quantities cannot vanish in general, neither 
can the pressure. Thus, the virial theorem which is combined 
with, but in no way related to the non-linear density evolution 
requires a pressure profile, even if the latter is based on cold 
dark matter. 



XII. CONCLUSION AND OUTLOOK 

We propose a way to set up a fully relativistic method 
to obtain the virial radius and the virial overdensity for 



Nevertheless, we have to admit that the assumption of an instan- 
taneously appearing pressure profile at virial redshift is a highly 
idealized concept. 



the EdS and ACDM cosmology. Within the assumption 
of an approximately time-like Killing vector field of the 
FLRW metric, static solutions for perfect fluid spheres in 
general relativity (namely the Tolman-Oppenheimer -Volkoff 
equation) have been successfully applied to extend the virial 
theorem in a consistent manner. The result has been inserted 
into the virialization equation which can be solved for the 
virial radius to find the corresponding virial overdensity. It 
turns out that the solution of Wang & Steinhardt (1998) 
([1]) for the virial radius in a ACDM cosmology is almost 
perfectly reproduced by our formalism which can also be 
shown analytically by performing the weak field limit. The 
first order post-Newtonian expansion has been investigated 
and the leading order corrections have been worked out 
and calculated numerically. We found out that they have a 
relative contribution of 10 _3 % with respect to the classical 
term which is of very small, but expected size. Although 
these corrections are of limited astrophysical interest, the 
concept itself is a small step towards a more fundamental 
understanding of virialization of spherical halos in the 
presence of DE. In addition, an iterative method has been 
set up to calculate the exact virial redshift numerically. The 
results of Lee & Ng (2010) ([6]) are reproduced extremely well. 

Naturally, galaxies and galaxy clusters are far more complex 
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than homogeneous and isotropic spheres, but the spherical 
collapse model provides a very simple semi-analytic method 
that already suffices to estimate important parameters like 
the virial radius and virial overdensity. The process of 
virialization itself is an additional condition that has been 
introduced to prevent a spherical overdensity of pressureless 
dark matter from collapsing into a singularity. A pressure 
profile, which a relaxed continuous spherical object must 
obey due to general relativity, can provide new insight into 
the process of virialization itself. 

There are certainly some topics this paper cannot address, 
because they are far beyond its scope. Although the spheri- 



cal collapse models are powerful tools to obtain estimates of 
the evolution of structures in the universe, they are limited 
by their simplicity. In particular, the fact that the spheri- 
cal overdensity is in no way embedded continuously into the 
background Friedmann universe is still very idealized and 
dissatisfying. Secondly, DE cannot be described yet in a 
self-consistent way with general relativity, since local energy- 
momentum-conservation, which is required by the theory, is 
not fulfilled and a coordinate representation of the two fluids 
is missing. Approaches based on the Lemaitre-Tolman-Bondi 
models (see [31-36]) as well as the presented work of Mis- 
ner and Sharp (1964) (see [26]) are promising candidates for 
following investigations. 
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Appendix A: Approximately time-like Killing vector 
field of FLRW spacetime on halo scales 

Let us start with the FLRW-spacetime given by the metric 



ds = -dr + 



(dr 2 + r 2 dO, 2 ) . (Al) 



This model is based on isotropy and homogeneity of the 
3-dimensional spacelike hypersurfaces describing a space of 
constant curvature k. 

Isotropy requires the existence of a coordinate frame in which 
spatial rotations are isometries. Given that frame, the most 
general ansatz for a Killing vector field of Eq. (Al) is: 



Let us now apply the following approximation: 

Since the spherical overdense region is of the scale of 1 — 10 
Mpc which is small compared the Hubble radius Rh ~ ~ 
1 — 3 Gpc, it can be assumed that 



2 1 

r «k 



which implies 



2 



1 + 



kr- 



< l. 



(A15) 



(A16) 



Thus, we can write approximately for Eq. (A13) (adding the 
factor c again) : 



K = A{r,t)d t + B(r,t)d r (A2) 

with A and B being arbitrary scalar functions of radius and 
time. 

Eq. (A2) has to be inserted into the Killing equation to obtain 
any relation between A and B. 

(£*3) M „ = 0. (A3) 
The Lie derivative of a rank (0, 2) tensor can be written as 

K x dxg^ + gxud^K x + g^d„K x = 0. (A4) 

Inserting Eq. (A2) into Eq. (A4) the following three con- 
straints can be obtained: 



(£k-s) 



2A = 



(A5) 



{C K g)^ ■ (2aA + 2aB') (l + ^r 2 ^ - kraB = (A6) 
{C K g) 22 : {2aAr 2 + 2arB) (l + \r' 



(£ K g) c 



kr 6 aB = 
A' 



(1 + T) 2 



(A7) 
(A8) 
(A9) 

(A10) 



(C K g) :i3 = sin 6(£ K g) 22 = 
Eq. (A5) implies that A has no time-dependence, thus 



A(r,t) = A(r). (All) 

If Eq. (A6) and Eq. (A7) are combined, we can obtain the 
differential equation in B 



-B 



= —B 



-B 



Hr 

c 
c 

Hor = 



-B 



E{a) 
Rh 



(A17) 



because E(a) = \[^ar 3 + + Q, { ° ] a~ 2 ~ O(l) (within 
the range of parameters we consider) and Rh = c/Hq. 



(A18) 



This implies 

B _ r 
~A~ R~h~' 

In our considered spherical collapse model, parameters are 
given within the following range 

o c ,o ta ~O(l)-O(10- 1 ) 
r ~ 1 - lOMpc 
Rh ~ 3Gpc. 



Inserting this, we end up with 

\B\ 



\A\ 



< 3.3 • 10" 



0.3% 



(A19) 



which means that the Killing vector field is approximately 
time-like: 



K « A(r)d t 



(A20) 



B' 



(A12) 



This can be inserted into Eq. (A6) to obtain a relation be- 
tween A and B 



A = 



B 
rH 



2 



fcr 2 



(A13) 



It has to mentioned for completeness that Eq. (A13) has 
to be consistent with Eqs. (A5) and (A9). Calculating the 
derivative of Eq. (A13) with respect to r and equating this 
with Eq. (A9) leads to an additional constraint for B: 



Appendix B: Derivation of the relativistic potential 
energy 

The following derivation is essentially taken from N. Strau- 
mann ([21]) and will be only slightly modified for our 
purposes: 

Let us consider a particle representation of dark matter and 
compare the total mass M of the spherical halo with the rest 
mass of the gravitationally interacting dark matter particles. 
The total rest mass of all DM particles is given by 



A' 



B 
rH 



kr 



a 2 B 



B 



kB 
a?H' 



M = Nmp 



(Bl) 



with being the total number of particles and ttln the rest 
(A14) mass of a single particle. Let us define a current density J as 
a one-form such that we can express the number of particles 
via a surface integral: 
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N - 



L 



* J. 



J can be expressed as 



The proper internal energy can be decomposed into a total 
(B2) kinetic and a total potential energy of the system such that 
T + V — M — M . Let us insert the integral for the particle 
number given by Eq. (B7) 



J= J a 



*J=J^ * 9" = J u rf with rf = *6 



(B3) 



with respect to an arbitrary dual basis {# M }. 
We transform the integral by evaluating the Hodge dual ex- 
plicitly: 21 



/ * J = / Jpff 

J t = const J t — const 

= / w 

J t = const 

= f Jo* 9". 

J t — const 



(B4) 



Assuming the metric ansatz given in Eq. (4.1) and defining 
the dual basis like 



9° = e a dt, 9 1 = e b dr, 9 2 = r 2 d9, 9 Z = r 2 sin 2 
*9° can be evaluated: 



(B5) 



f R 

m N N — / 4irr 2 e b m N n(r)dr 
Jo 

f R 

— / (p(r) — e(r)) 4ttt e dr 
Jo 



(B12) 



M -T-V 

r-R 



/ 4nr 2 p(r)dr -T-V 
Jo 



where we have assumed that the total mass is simply the vol- 
ume integral of the density profile 



f R 2 
M — Anr p(r)dr. 

Jo 

Solving Eq. (B12) for T + V, we obtain 

T+V = J 4-Kr 2 e b e(r)dr + J 4nr 2 p(r)( 
This leads to the definition 



(B13) 



1 - e 



')dr. (B14) 



*9° = e a * dt 



er 2 sin 9dr A d9 A • 



(B6) 



Since, in a static configuration, the 6" 1 are an orthonormal 
system, the above result can be expected. 

Thus, we are left with 



N 



f J06I 1 A 9 2 A 9 3 

J t=const 

L 
L 



Joe b r 2 sin 9dr A d9 A < 



t — const 
R 



(B7) 



4nr 2 Joe b dr. 



The number density n(r) can be obtained by projection of 
J onto the four velocity tt M being (1,0,0,0) T in the chosen 
coordinate frame: 



n(r) = -u^J^ = J 



such that we get 



N 



4nr 2 n(r)e b ^ dr 



(B8) 



(B9) 



We can define the proper DM energy density (total energy 
density with subtracted particle rest energy density) 



4nr 2 e(r)e b dr 



Jo 

U = j 4nr 2 p(r) ^1 — e b ^ dr. 



(B15) 
(B16) 



Consider a top-hat density profile and a two component fluid 
consisting of DM and a cosmological constant such that 



e 2b = 



r -i^ with A= *^l( p + PA ). (B17) 



Thus, we finally end up with: 



A 2 ' 

4nr £ 



-r 

Jo 

f R 2 ( 

= / 4nr p I 1 



VI - Ar 2 



dr 



Vl - Ar 2 



dr. 



(B18) 
(B19) 



Appendix C: Static, spherically-symmetric field 
equations with homogeneous DE 



Consider a two component fluid described by 



= (P+P) U H U V + P9nv 



TjS ] = {pq +pq)u u u v +p Q g llv 



(CI) 
(C2) 



e(r) = p(r) - m N n{r) (BIO) 
which corresponds to an intuitive proper internal energy of 



E = M-M = M- iVmjv 



(Bll) 



T — T< m ) i t(Q) 



(C3) 



where the densities p and pq are assumed to be constant 
and the quintessence component has an equation of state 
PQ = w>pQ with constant w. Energy-momentum conservation 
is fulfilled separately for each fluid component: 



21 Since we consider a static configuration we can find a coordinate 
frame in which the J components vanish. We will assume that 
in the following without loss of generality. 



V»T» v ' m = 
v T nv,Q = Q 



(C4) 
(C5) 
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The static, spherically symmetric field equations of this set-up 
are 



Guv = 8tyGTuv 



1 2b' 



= 87rG(p + p Q ) (C6) 



a" -a'b' + a' 2 + a ^ - 



= 8tvG (p + wp Q ).(C8) 

We have to find out whether there are conditions for the 
solvability of the field equations without using any concrete 
solution for a and 6 such that effects resulting from boundary 
conditions are excluded. 

With the help of Eq. (C6) we can express b 1 : 

b' = y r (l-(l-8nG(p + p Q )r 2 )e 2b ). (C9) 
In the same way Eq. (C7) can be solved for a': 

a = - ^ (l - (1 + 8ttG (p + wp Q ) r) e 26 ) . (CIO) 
If we add Eq. (C6) and Eq. (C7), we will get 



a = —b' + 4nG (p + p + pq(l + w)) re 2 



(Cll) 



If Eq. (C9) is inserted into that expression, we will obtain 
Eq. (CIO) so the first two field equations are consistent. 

Energy-momentum conservation for the matter component of 
the fluid means 

V M T Ml/ ' m = 0. (C12) 

Projecting this onto the space perpendicular to the velocity 
flow, we obtain the relativistic Euler equation 



{g av + u a u v ) V^T"" =0 



which becomes 



(p + p) V u u = — gradp — uV u p. 
In case of a static configuration, we are left with 



(C13) 
(C14) 



-p' = a'(ft + p) (C15) 

which is basically the hydrostatic equilibrium condition for 
the matter configuration of our system. 

Consider the derivative of Eq. (CIO) 



a" (l + e 2b {-1 + 8ttG (p + wp Q ) r 2 + 8irGp'r z + 2b' r 

{l + 8nG(p + wp Q ))r 2 }) . 

(C16) 

Inserting Eq. (C9), Eq. (CIO) and Eq. (C15) into Eq. (C16) 
and simplifying leads to 



-^e 4b (l + 8nG(p + wp Q )r 2 ) 
(1 - 4ttG (p - p + 2 PQ ) r 2 ) . 



(C17) 



On the other hand a" can be expressed with the help of Eq. 
(C8) 



/,/ 12 a 

a b — a 



1/ 



+ 8tyG (p + wp Q ) e 



2b 



(C18) 



If we plug in Eqs. (C9) and (CIO), we can obtain (after some 
algebra) 



a " =2^ + 2^ e2i> t 47rG (P + PQ) r " + 207VG {P + WpQ) f ' 2 ] 
1 e 4b (1 + 8nG (p + wp Q ) r 2 ) (1 - 4ttG (p + p Q 



2r 2 
p - wp Q )) 



(C19) 



For Eqs. (C6) - (C8) being consistent, Eq. (C17) and Eq. 
(C19) have to be identical which means that each coefficient 
belonging to the same order of e 2b has to be equal for each 
radius r. 



• 0th order: 



1 1 

272 ~ 



trivially fulfilled 



1st order: 



^-47rG (5p + 4wp Q + p) r 2 



~2r 2 



[4nG (p + pg) r 2 + 20ttG [p + wpq) r 2 ] 

5p + p + 4wpQ = p + pq + 5p + 5wpQ 
=> Aw — (1 + 5w) 
=> w — -1 



2nd order: 
1 



2r 2 
1 

= 2r2 



(1 + 8ttG (p + wp Q ) r 2 ) (1 - 4ttG (p - p + 2p Q ) r 2 ) 
(l + 87rG (p + wpq) r 2 ) (l — 47rG (p — p + pq - wpq) r 2 



If we require that p =fc —wpQ — l/(87rGr 2 ) which we have to 
assume in the general case, we can say 



4ttG (p - p + 2p Q ) = 4ttG (p - p + , 

2p Q = p Q (1 - W) 

w = -1 



wpq) 



Thus, Eq. (C6) - Eq. (C8) necessarily require w 
order to be consistently solvable. 



-1 in 



Appendix D: Iterative method to find the virial 
redshift 

Consider Eq. (2.1) and solve this for the non-linear density 
contrast 5(a). Once this is done, Eq. (2.2) can be used 

A(r,o) = ^=cfrV = 1 + (Dl) 



(D2) 



Pb(a) \y 
This allows to express y(a) with the help of S 

1/3 



v(a) 



c 



fflta \l + <5(a) / 

Using Eq. (D2) and the virialization equation (Eqs. (2.3), 
(7.1), (8.23)), an iterative method can be constructed to find 
the virial scale factor a v ; r . 
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Starting from a' - 1 = a c , we will proceed in the following way 



(0) ( 2 -3) (0) (D2) 

(1) (2^ (!) (£2) 

(2) (2.3) (2) (D2) 

1 —> y J — 



« (l&>) 

« (&>) 



,(2) 



(D3) 



r' / (")\ 



(n) (2^3) (n) (£21 



a 



(n + 1) 



In each step, the quantity a(y^) is needed which is given by 
the root of the equation 



V(a) -»«=()—► a(y<g) 
and which can be found numerically in each step i. 



(D4) 



It turns out that this method converges 22 , so that the follow- 
ing condition can be applied to stop the iteration: 



("+ 1 ) _ a (") 



(D5) 



with an appropriate tolerance ioL Thus, the final result can 
be defined as: 



- 1. 



(D6) 



It will turn out that the virial overdensity changes significantly 
if the exact quantities are evaluated at z = z v i r : 
Let us calculate Ay at redshift z c = in the EdS model 
which corresponds to a virial redshift of z v -„ = 0.065 and a 
turn-around redshift of z ta = 0.587. 

• Case 1 (z — z c ): 

A V =(?)VW(^) 3 = 177.518 (D7) 

• Case 2 (z — 2 v ir): 



1 + Zta \ 3 / J_\ 3 = 
1 + z vil J \ y vil J 

This has also been predicted analytically by Lee, Ng (2010) 
(see [6]). 



22 Given a tolerance of tol ~ 10 7 , the iteration can be stopped 
within 3 to 4 steps. 



